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Abstract 

In this paper we study the su(to) spin Sutherland (trigonometric) model of Dn type and 
its related spin chain of Haldane-Shastry type obtained by means of Polychronakos's 
freezing trick. As in the rational case recently studied by the authors, we show that 
these are new models, whose properties cannot be simply deduced from those of their 
well-known BCn counterparts by taking a suitable limit. We identify the Weyl- invariant 
extended configuration space of the spin dynamical model, which turns out to be the N- 
dimensional generalization of a rhombic dodecahedron. This is in fact one of the reasons 
underlying the greater complexity of the models studied in this paper in comparison with 
both their rational and BCn counterparts. By constructing a non-orthogonal basis of the 
Hilbert space of the spin dynamical model on which its Hamiltonian acts triangularly, we 
compute its spectrum in closed form. Using this result and applying the freezing trick, we 
derive an exact expression for the partition function of the associated Haldane-Shastry 
spin chain of Djv type. 
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1. Introduction 

Recent studies have revealed that exactly solvable and integrable one-dimensional 
quantum many body systems with long-range interactions [lHSl are closely connected 
with a wide range of topics in modern physics as well as mathematics. In particular, 
this type of exactly solvable systems have appeared as prototype models of various con- 
densed matter systems exhibiting generalized exclusion statistics [cj-lioj]. quantum Hall 
effect 11] and quantum electric transport phenomena 12, 13| ■ In the context of high- 



energy physics, the dynamics of particles or fields in the near-horizon region of black holes 
has been described through such integrable systems More recently, quantum in- 

tegrable spin chains with long-range interaction have played a key role in calculating 
higher loop effects in the spectra of trace operators of planar N = 4 super Yang-Mills 
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theory 17-13]. Furthermore, this type of quantum integrable systems are found to be 
connected with different areas of mathematics like random matrix theory [20j, multi- 
variate orthogonal polynomials 21 - 24 1 , Dunkl operators [2|| [26[ , and Yangian quantum 
groups 



27H30(. 

Due to such a large variety of potential applications, the construction of new quantum 
integrable systems with long-range interactions and the computation of their exact solu- 
tions have emerged as an important area of activity in the current literature. The study 
of this type of systems with dynamical degrees of freedom was pioneered by Calogero, 
who found the exact spectrum of an iV-particle quantum system on the line with two- 
body interactions inversely proportional to the square of the distances and subject to a 
confining harmonic potential 1]. An exactly solvable trigonometric variant of this ratio- 
nal Calogero model was subsequently proposed by Sutherland 0, H| ■ The particles in this 
so-called Sutherland model move on a circle, with two-body interactions proportional to 
the inverse square of their chord distances. In a parallel development, Haldane and Shas- 
try found an exactly solvable quantum spin-i chain with long-range interactions [H, Q. 
The lattice sites of this su(2) Haldane-Shastry (HS) spin chain are equally spaced on 
a circle, all spins interacting with one another through pairwise exchange interactions 
inversely proportional to the square of their chord distances. A close relation between 
the HS chain with su(m) s pin degrees of freedom and the su(to) spin generalization of the 
Sutherland model [3ll - t33| was subsequently established by using the so-called "freezing 
trick" 0, Hil ■ More precisely, it is found that in the strong coupling limit the particles in 
the spin Sutherland model "freeze" at the coordinates of the equilibrium position of the 
scalar part of the potential, and the dynamical and spin degrees of freedom decouple. 
The equilibrium coordinates coincide with the equally spaced lattice points of the HS 
spin chain, so that the decoupled spin degrees of freedom are governed by the Hamilto- 
nian of the su(m) HS model. Moreover, in this freezing limit the conserved quantities 
of the spin Sutherland model immediately yield those of the HS spin chain, thereby 
explaining its complete integrability. Application of this freezing trick to the rational 
Calogero model with spin degrees of freedom leads to a new integrable spin chain with 
long-range interaction [7J . The sites of this chain — commonly known in the literature as 
the Polychronakos or Polychronakos-Frahm (PF) spin chain — are unequally spaced on 
a line, and in fact coincide with the zeros of the Hermite polynomial of degree N [35| . By 
applying the freezing trick, the exact partition functions of the PF and HS spin chains 
have also been exactly computed [3H l37| . 

The above mentioned type of quantum integrable systems can be generalized to form 
a much wider class by taking advantage of their hidden mathematical structure. Indeed, 
Olshanetsky and Perelomov established the existence of an underlying Am root system 
structure for both the spinless Calogero and Sutherland models, and constructed gener- 
alizations thereof associated with any (extended) root system [4| . Spin generalizations of 
the BCn Calogero and Sutherland models have also been proposed, and various prop- 
erties of the related lattice models of HS type have been studied with the help of the 
freezing trick 3^-43|. Among the other classical root systems, the exceptional ones are 
comparatively less interesting, since their associated models consist of at most 8 parti- 
cles. Until recently, the Bn, Cn and Dm Calogero-Sutherland models (particularly the 
corresponding spin models) have been largely ignored, probably due to the fact that they 
were believed to be simple limiting cases of their BCn counterparts. However, in a recent 
paper [44j the present authors have computed the spectrum of the su(m) spin Calogero 
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model of Dm type, thereby showing that this model is in fact a singular limit of its BCm 
version. More precisely, it is well known that the Hilbert space of the spin Calogero 
model associated with the BCm root system can be constructed from the Hilbert space 
of an auxiliary differential-difference operator by using a single projector. In contrast, 
it is found that two independent projectors of BCm type with opposite "chiralities" are 
needed to construct the Hilbert space of the .DAT-type spin Calogero model from that of 
the corresponding auxiliary operator [44[ . Consequently, the Hilbert space of the latter 
model can be expressed as a direct sum of the Hilbert spaces associated with two differ- 
ent BCm models with opposite chiralities. This explains why the spectrum of the Dm 
model cannot be obtained as the limit of its BCm counterpart when one of the coupling 
constants tends to zero. In Ref. [44| we also studied the spin chain associated with the 
-Djv-type spin Calogero model, showing that its Hamiltonian differs from the limit of its 
BCm analog by a term which can be interpreted as an impurity interaction at one end 
of the chain. By applying Polychronakos's freezing trick we were also able to compute 
the chain's partition function in closed form, showing again that its spectrum markedly 
differs from that of its BCm counterpart. 

In this paper we study the trigonometric variant of the Djy-type spin Calogero model 
and its freezing limit, i.e., the Djv-type spin Sutherland model and its related spin 
chain. Just as in the rational case, the Hamiltonian of the Dj\r-type spin Sutherland 
model can be formally obtained as a certain limit of its BCm counterpart. However, the 
relation between these models turns out to be even more subtle than in the rational case. 
Roughly speaking, this is due to the fact that the Weyl-invariant extended configuration 
space of the Dm model — which turns out to be the iV-dimensional generalization of a 
rhombic dodecahedron — does not coincide with that of the BCm model, which is simply 
a hypercube. As a consequence, the (scaled) Fourier basis of the Hilbert space of the 
BCm model's auxiliary operator no longer spans a complete set of the Hilbert space of the 
auxiliary operator of the Dm model. This entails an additional level of difficulty (but also 
of interest) by comparison with the rational case, for which the auxiliary operators of the 
BCm and Dm models share the same Hilbert space. On the other hand, as in the rational 
case, we shall still need two projectors of BCm type with opposite chiralities in order to 
construct the Hilbert space of the Dm spin model from that of its corresponding auxiliary 
operator. Therefore, the Hilbert space of the Dm spin model actually consists of four 
— and not two, as in the rational case — different sectors, characterized by their chirality 
and parity under reflections of the particles' coordinates. This fundamental difference 
explains why the spectrum of the Dj\r-type spin Sutherland model is essentially different 
from that of its BCm counterpart. It also accounts for the greater complexity of the 
partition function of the associated chain of Dm type (which we have also computed 
in closed form by means of the freezing trick) compared to its BCm version studied in 
Ref. Q. 

The paper is organized as follows. In Section [5] we introduce the Hamiltonians H, 
H sc and % of the D/v-type spin Sutherland model, its scalar version, and the associated 
spin chain of HS type, respectively. We show that the sites of this chain, defined as the 
coordinates of the (unique) equilibrium point of the scalar part of the spin Hamiltonian 
in the principal Weyl alcove of the Dm root system, can be expressed in terms of the 
roots of a suitable Jacobi polynomial. Using this characterization, we prove that the 
Hamiltonian % differs from the limit of its BCm counterpart by a spin reversing term at 
each end of the chain. 
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Section [3] is devoted to the computation of the spectrum of the Hamiltonians H and 
H sc using an auxiliary scalar differential-difference operator H' . In order to improve 
the clarity of the exposition, we have divided it into four subsections. In the first one, 
we show that for any one-dimensional representation tt of the Dn Weyl group 2B, the 
Hamiltonians H and _ff sc are equivalent (isospectral) to their 7r-symmetric extensions 
to the 2U-orbit C of the configuration space. The representation it is then uniquely 
determined by requiring that the action of the extension of H (resp. H sc ) coincide with 
that of H' (g) II (resp. H') on the subspace of L 2 (C) ® £ (resp. L 2 (C)) of 7r-symmetric 
functions, where £ denotes the su(m) spin space. This property, together with the fact 
that H' commutes with the projector onto 7r-symmetric functions, enable us to evaluate 
the spectra of H and H sc by triangularizing the simpler operator H'. The first step 
in this direction is to construct a (non-orthogonal) basis of L 2 (C), which includes as a 
proper subset the limit of the basis of L 2 ([— ^, ^]) for the BCn model. As a by-product 
of this construction, we show that the translations of C generate a tessellation of the 
TV-dimensional Euclidean space. By expressing W as a sum of squares of a suitable 
family of Dunkl operators, we show that when the above basis of L 2 (C) is ordered in 
an appropriate way the action of H' becomes triangular. Using this basis, in the last 
subsection we construct a corresponding (non-orthogonal) basis of the subspace of ir- 
symmetric functions L 2 (C) (8 S (resp. L 2 (C)) in which the action of H (resp. H sc ) is 
also triangular. This completes the calculation of the spectra of the Hamiltonians H and 

H sc . 

In Section 2] we compute the partition function Z of the spin chain H using Poly- 
chronakos's freezing trick. More precisely, we evaluate the partition function Z as the 
large coupling constant limit of the quotient of the partition functions of H and _ff sc . 
The resulting formula exhibits a greater complexity than its rational counterpart and, 
in particular, cannot be expressed in a simple way in terms of the partition functions of 
the BCn trigonometric spin chains. We end up this section by showing that the latter 
expression for the partition function is indeed a polynomial in q = Q~ 1 /( k ^ T ) ^ as should be 
the case for a finite system with nonnegative integer energies. Finally, a brief summary 
of the paper's results is presented in Section [5] 

2. The models 

Our starting point will be a brief review of the su(m) spin Sutherland model of BCn 
type, with Hamiltonian [39|, |4^| 

H (B) = - E d l + a E [ sin ~ 2 x ij (« + s ^ + sin ~ 2 4 (° + § v)] 

+ b nm~ 2 Xi {b - eSi) + 6' E cos " 2 ^ ( b> ~ eS i) ■ i 1 ) 

i i 

Here the sums run from 1 to N (as always hereafter, unless otherwise stated), a,b,b' > 
1/2, e = ±1, xfj = Xi ± Xj. The operators Sij and Si in the latter equation act on the 
finite-dimensional Hilbert space 

S = ^ |si, ...,s N ) \si = -M, -M + 1, . . . , Mj, M = — ^— €-, (2) 
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associated to the particles' internal degrees of freedom, as follows: 

Sij \s\ , . . . , Si, . . . , Sj, . . . , Sjv) = \si, . . . , sj , . . . , Si, . . . , sjv) j .g. 
Si\si, . . . , Si, . . . , s N ) = |si, . . . , -Si, . . . , s N ) . 

We have also used the customary notation Sij = SiSjSij. Note that the spin operators 
and Si can be expressed in terms of the fundamental su(to) spin generators Jg at 
the site k (with the normalization tr(J^J^) = \8 ai ) as 

1 m 2 -l 

Sij = - + 2 V J?Jf , Si = y/tol J\ . 

a — 1 

Due to the singularities at the hyperplanes X4 ± Xj = kit, Xi = kir and Xj = \ + (with 
1 ^ i < J ^ N and fc 6 Z), the configuration space of the Sutherland Hamiltonian (JXJ) 
can be taken as the principal Weyl alcove 

A {B) = jx e R N : < xi < x 2 < ■ ■ ■ < x N < || (4) 

of the BCn root system. The spectrum of the spin model ([1} was computed in Ref. [42[ 
by constructing a suitable basis of its Hilbert space in which the Hamiltonian is 
represented by a triangular matrix. 

Applying the so-called freezing trick [7J to the Hamiltonian (Q} with b — /3a and 
b' = j3'a one obtains the su(m) Haldane-Shastry (antiferromagnetic) spin chain of BCn 
type, whose Hamiltonian we shall take as 

+ \ H (p sin ~ 2 e * + p' cosr2 O ( x - eS ^ ■ ( 5 ) 

Here Ofj = 9i ± where = (#1, . . . , #at) is the unique equilibrium [4l| in the set 
of the scalar potential 

c/ (B) (x) = J2 ( sin_2 + sin ~ 2 x ij) + (^ 2 sin ~ 2 Xi + p' 2 cos ~ 2 x ■ ( 6 ) 

As shown in the latter reference, the lattice sites Oi are related to the zeros Q of the 
Jacobi polynomial pj^ ^ by 

Q = coB(2fl() . (7) 

The spin chain ([3]) was studied in Ref. (42[, where its partition function was computed 
in closed form with the help of the freezing trick. 

The Hamiltonian H of the su(m) spin Sutherland model of D n type is defined by 
setting b = b' = in Eq. (JTJ, i.e., 

H = - E d l + a E [ sin_2 (° + s ^ + sin ~ 2 4 ( a + 5 « )] ■ ( 8 ) 

i i=jLj 
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The configuration space A of the D n model (O is determined by the hard-core singular- 
ities of the Hamiltonian on the hyperplanes Xi ± Xj = kir (with i ^ j and k € Z). More 
precisely, we shall take as A the open subset of R N defined by the inequalities 

< Xi ± xj < 7r , 1 < j < i ^ N. (9) 

If N > 2 (as we shall assume hereafter), it is straightforward to check that this set can 
be equivalently expressed as 

A = {X E R N : \xi\ < X 2 < ■ ■ ■ < X N < 7T - x N -i} , (10) 

which is again the principal Weyl alcove of the Djy root system 

-(±e i ±e j ), l^i<j^N. (11) 

The points x £ A clearly satisfy the inequalities 

< x 2 < ■ ■ ■ < x N -i < 7r/2, xi > -7r/2, x n <tt; (12) 

note, in particular, that the set A properly contains the BCn configuration space 

Similarly (cf. (J5J) and ©), wc define the Hamiltonian of the su(m) HS spin chain of 
Dn type as 

n = \ E [ sin_2 *V + + sin_2 *ii 0- + &4 > (13) 

where the lattice sites i9j are the coordinates of the unique minimum i? in the set A of 
the scalar potential 

(7(x)^(sin- 2 xr.+sin- 2 x+). (14) 

Heuristically, the relation between the spin dynamical model (O and its associated spin 
chain (|13|) can be explained as follows. Defining the coordinate-dependent matrix multi- 
plication operator 

i<j 

the spin Hamiltonian ((SJ) can be decomposed as 

F = ff sc + 4a/i(x), (15) 



where 



H sc = -Y,dl+a(a-l)U( X ) (16) 



is the Hamiltonian of the scalar Sutherland model oi Dn type. Thus, for sufficiently 
large a all the eigenfunctions of H sc are sharply peaked around the unique minimum # 
of the scalar potential U in the set A [iBj]. Hence, if ^i(x) is an eigenfunction of H sc 
with energy Ef c and |oj-) is an eigenstate of the chain H with eigenvalue £j, for a ^> 1 
we have 

/l(x)<Pi(x)|(Tj) ~ <£j(x)/l(#)|<7j) = ipi(x)H\<7j) = £j<Pi(pt)\(Tj) . 
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By Eq. (|15|) . H is approximately diagonal in the basis with elements ^(x)|<7j), and its 
eigenvalues Eij satisfy 

E tj ~ E s t c + 4a£j , a > 1 . (17) 

It was shown in Ref. [4l| that the scalar potential U (x) has a unique minimum in the 
configuration space A, which coincides with the unique maximum in this set of the ground 
state wave function of the scalar Sutherland Hamiltonian of Dn type (|16[) . given by 



/o(x) = J^J | sin i 4J sin | a . (18) 

i<j 



The lattice sites i9< of the chain ([13)) are thus the unique solution in A of the nonlinear 
system 

(cottfr. + cottft.) =0, lsCisCA^. (19) 

As in the -BCjv case, we define the variables £j by 

& = cos(2i? 4 ) , 1 < i < iV . 
Note that, since i9 € A, we obviously have 

1 > 6 >&>••• > 6v-i > &v > -1 ■ (20) 
In terms of the new coordinates £j, the system (fl9l) can be written as 



t 1 "^) E F3T = °> Ki<N. (21) 



Since £i — £j > for all j > 1 and £jv — £j < for all j < AT, from Eq. (1211) it immediately 
follows that £f = ^ = 1, so that £i = -£at = 1 by Eq. ([20]). Substituting into ([HJ we 
obtain the following system for the remaining coordinates £2, ■ • ■ , £n—i- 

AT-l 1 

( 1 -^E7T7 =2 ^' 2<i<JV-l. (22) 
3=2 5i ^ 

Note that the latter system is invariant under the transformation £j 1— >• — £j, so that (by 
uniqueness) £j = Cw+i-i- Comparing (|22[) with the system 

N> 

^-^Et-Z^^' + C + W- (23) 

satisfied by the zeros Ci (i = 1, • • • , N') of the Jacobi polynomial Pj^, (cf. Ref. (iij). 

we conclude that the coordinates £2, ■ • ■ , 6v-i are the zeros of P^l- (Note that Pyl_2 is 
proportional to the Gegenbauer polynomial cf. Ref. 0.) In terms of the original 

site coordinates §i we have 

= 1? l <^ 2 <---<^ iV _l<^ iV =|, 
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with P^l ( cos(2i?i)) = for % = 2, . . . , N - 1. Note that Eqs. ([2"T ]l -([2"3 |) also yield an 
alternative characterization of the coordinates £j as the N roots of the Jacobi polynomial 
P N l ' , which was to be expected, since the potential C/( B ' in Eq. (jSJ) reduces to the 
Dn potential U when /3 = j3' = 0. The equivalence of both characterizations of the 
site coordinates is easily established with the help of the identity AP^' (t) = {t 2 — 
l)P^l(t),d. @. 

We shall next discuss the precise relation between the Dm spin chain Hamiltonian (fT5|) 
and the limit as (f3, (3') — > of its BCm counterpart ([5]). To this end, we use the 
trigonometric identities 

sin~ 2 6i = - - - , cos -2 



1-C; 1 + C»' 

and note that as (/3,/3') — > all the roots £, of the Jacobi polynomial P^f -1 '^ tend 
to the corresponding roots £j of Pjy '~ . Thus all the terms in the last sum in the 
Hamiltonian (O tend to zero as (j3,j3') — > (0,0), except the first and the last one. In 
order to evaluate the limit of these two terms, we divide (j2"3")l by 1 ± Q and sum the 
resulting equation over i, obtaining 

2/3 E rh~ = 2/3 ' E ttt- = n(p+p+n-i). 

i i 

Hence 

lim = lim 2/3 = 7V(iV - 1) . (24) 

(/s,/J')->-° 1 - Ci WHO 1 + Cw 

Since £i & as {P, P') and 0,-, -di & A (recall that C A), we have lim^^^g 6>j = 
&i for all i = 1, . . . , N. From Eqs. ©, (O and flU) it immediately follows that 



1 



(/3,/3')~>0 2 



lim H (B) = H + ~ N{N - 1) 1 - + ^w) . (25) 
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Thus the limit as (/3, /3') — > of the Hamiltonian of the HS chain of PCV type yields its 
Dn analog, plus an additional term which can be interpreted as an "impurity" at both 
ends of the latter chain. 



3. Spectrum of the dynamical models 

The aim of this section is to compute the spectra of the su(m) spin Sutherland model 
of Dn type ([SJ and its scalar counterpart ([T5| . In order to facilitate this computation, 
we introduce the auxiliary scalar operator 

H ' = - E d l + a E [ sin ~ 2 x ij (° - K « ) + sin_2 4 ( a - KiJ )] . ( 26 ) 

i ijtj 

where K+j and Ki are coordinate permutation and sign reversing operators, defined by 

{Kij f)(xi, . . . , Xi, . . . , Xj, . . . , Xjv) = f{x\, . . . , Xj, . . . , Xi, . . . , xn) , 

(Kif)(xi,...,Xi,,..,x N ) = f(xi,...,-Xi,...,x N ), 
and JQj = KiKjKij. 



3.1. Extensions of H and H sc 

Due to the character of their singularities, the operators H and H sc are naturally 
defined on suitable dense subspaces of the Hilbert spaces L 2 (A) <g) E and L 2 (A), respec- 
tively. On the other hand, the appearance in the RHS of Eq. (OBI) of the generators KiKj 
and Kij of the Weyl group 2U of the Dm root system entails that the auxiliary operator 
H' is defined instead on a dense subspace of L 2 (C), where C = W(A). One of the key 
ingredients of the method we shall use consists in replacing the operators H and -ff sc by 
suitable equivalent (isospectral) extensions H and H sc thereof to appropriate subspaces 
of L 2 {C) (g> E and L 2 (C), such that H = H' ® I and £T SC = H' in the latter subspaces. 

We shall start by showing that the set C is explicitly given by 

C = {x £ R N : < \ Xi ±xj\ < 7T, f <i < j < N} , (27) 

a characterization which will prove useful in what follows. Recall, to this end, that 
2B is generated by coordinate permutations and sign reversals of an even number of 
coordinates [4!|. Since A is defined by the inequalities ((HJ) , it is obvious that C C C*, 
where C* denotes the RHS of Eq. ([2"7]). Thus, we need only show that C* C C, or 
equivalently, that for every x* £ C* there is an element W £ 2U such that Wx* £ A. To 
this end, we shall make repeated use of the following elementary fact: 

xeC*, \xi\ < x 2 < ■ ■ ■ < x N =^ xeA. (28) 

Suppose, then, that x* £ C* . By reversing the sign of an even number of coordinates 
and applying a suitable permutation, we can transform x* into another element y £ 
W{C*) = C* satisfying 

Vi<V2<---<Vn and > , i = 2,...,N. 

If 2/1+2/2 > 0, then \yx\ < yi < ■ ■ ■ < yn, and hence y £ A by Eq. (|2"51) . Suppose, on the 
other hand, that y\ + j/2 < 0, so that 

yi < < 2/2 < • • • < yN ■ 

Calling z\ — — 2/2> z 2 = — 2/1 and Zi = yi for i ^ 3, we have z £ C* and 

zi +z 2 = -(2/1 + 2/2) > 0, zi < z 2 |zi|<-?2- 

If Z2 < Z3, then z€ Aby Eq. (|2^|) . Otherwise, the inequalities z\ + Z3 = 2/3 — 2/2 > and 
z 3 - zi = 2/3 + 2/2 > imply that 

|zi| < z 3 < • • • < z 2 < ■ ■ ■ < z N . 

Applying a suitable permutation to z, we obtain a new element u £ C* satisfying |ui| < 
1*2 < • • • < un, which belongs to A again by Eq. 



Remark 1. Note that the analogous set for the BCn root system, which is simply 
the hypercube (— f , f ) N minus the singular hyperplanes Xi±Xj — 0, l^i^j^N 



(cf. Ref. [42(1). is clearly contained in C by Eq. (|27)) . 



Remark 2. For N = 3, the set C is a rhombic dodecahedron (a zonohedron with 12 
equal rhombic faces) centered at the origin [50|, with edge length \/3n/2. 
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Our next aim is to replace the operator H by an isospectral extension H thereof acting 
on suitably (anti)symmetrized wave functions defined on the set C. If the extension is 
appropriately chosen, we shall see that H = H' ® 1 on the Hilbert space of H, so that the 
computation of the spectrum of H reduces to the analogous (but considerably simpler, 
in practice) task for H' . 

Before proceeding with our construction, we need to introduce some additional no- 
tation. Given an element W of the Weyl group of Dm type 2U and a factorized spin 
function <p\s) G L 2 (C) <g> £, where |s) = \si, . . . , sn) is an element of the canonical spin 
basis, we define the action of W on ip\s) in the usual way: 

W(<p\s)) = (<po W-^lWs) . (29) 

Extending this definition by linearity to the whole Hilbert space L 2 (C)(S)E, we obtain an 
action of 2XJ (in fact, a representation) on the latter space. Let now it : 2U — > C denote 
a one-dimensional representation of 2U; of course, since 2D is generated by reflections, 
we have 7r(2D) C { — 1, 1}- The symmetrizer associated with 7r is the linear operator 
defined on L 2 (C) <g) £ by 

^ = m J E>* {w)Wi (30) 

1 1 wem 

where |9U| = 2 1 iV! is the order of 2D. By construction, we have 

$ G A 7r (L 2 (C) ® £), WeW W$ = 7r(W)$, (31) 

so that A,r is the projector onto states with well-defined parity tt(W) with respect to any 
transformation W in the Weyl group 2D. 

Given a factorized spin function \f r = ?/j|s) 6 ^ 2 (^4) 55 E, we define its it-symmetric 
extension ^> G A X (L 2 (C) ® E) by 

*(x) = tt(VK x ) VG^x) I W x s) , x G C , (32) 

where W x denotes the unique element of 2D such that W~ x € A. As usual, the action of 
is extended to L 2 (A) ® E by linearity. It is easy to see that \P is the unique extension 
of $ to C which has well defined parity ir(W) under any transformation W G 2D. In 
view of Eqs. ([2U)) - (|c>0")) . with a slight abuse of notation we can write 

* = |2D| • K(*Xa) , (33) 

where \A is the characteristic function of A. 

The extension ~ : L 2 (A) <£> E — > (L 2 (C) ® E) is an invertible linear operator, 
its inverse being the restriction operator ~ : K^{L 2 {C) ® E) —> L 2 (A) ® E defined by 
$ = Indeed, the linearity of the operator ~ is obvious. As to its invertibility, note 

first of all that if * G L 2 (A) ® E by Eq. d32J) we have § = since W x is the identity 
when x belongs to A. On the other hand, if $ G A 7r (L 2 (C) <g> E) using Eqs. (JHOJ) and (jgg)) 
we obtain 

i=J] 7r(I¥)^($ X A)= J! tt(W)W($Xa)= J! ^W^)^"^ 1 ), 
wew wew wew 
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and hence, by Eq. (I3T1) . 



$(x) = $(x) Y, XA(^- 1 x) = <f>(x) Xyl (W^ 1 x) = $(x). 

wew 

Given a one-dimensional representation tt of 20 and a linear operator T acting on 
L 2 (A) ® £, it is natural to define its it -symmetric extension to the Hilbert space 
A,r(l/ 2 (C) ® S) by the prescription 

T,r $ = (T*)~ , $ G A T (L 2 (C) ® S) . (34) 

By the invertibility of the ~ operator, we have 

f OT = ~oTo(~)- 1 , 

so that the operators T and are isospectral. 

We now seek to find a suitable one-dimensional representation tt of 20 such that the 
7r-symmetric extension of H coincides with the restriction of H' ® II to A T (L 2 (C) ® S). 
In view of Eqs. (|8|)-(|26p. it suffices that tt satisfy 

KiKjh-T! = SiSjA-Tt , KijA-x = SijA.Tr . 

Hence 7r(T^) should be defined as the sign of the permutation part of W € 2U. From 
now on, when dealing with the spin Hamiltonian we shall take tt as above, and drop 
the subscript 7r from A ff and 

Turning now to the scalar Hamiltonian H sc , by Eqs. (|14p. (|16p. and (|26l) its extension 
i/ sc to the space A T (L 2 (C)) will coincide with the restriction of H' to the latter space 
provided that 

KiKjA n = K , KijK = K ■ 

Thus in this case tt(W) = 1 for all W G 20, so that for the scalar model (|lt>JI A T = A sc is 
the total symmetrizer with respect to both coordinate permutations and sign reversals 
of an even number of coordinates. 

By the previous discussion, we have reduced the problem of evaluating the spectrum 
of the Hamiltonians H sc and H to the analogous problem for the restrictions of the 
auxiliary operators H' and H' <g) II to the Hilbert spaces A SC (L 2 (C)) and A(L 2 (C) ® S), 
respectively. We shall next prove a more explicit characterization of these spaces that will 
be needed in the sequel. To this end, let be the projector onto states antisymmetric 
under particle permutations and with parity ±1 under reversals of coordinates and spins. 
We shall show that 

A(L 2 (C)®£) =A+(L 2 (C*)®£)®A-( J L 2 (C)<g)£). (35) 

Indeed, let A a denote the antisymmetrizer under particle permutations, and let {W i } 
be the set of reversals of an even (+) or an odd (— ) number of coordinates and spins. 
We then have 

^ i=l i=l ' 
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and hence 



2 



which establishes ((35)) . Similarly, if is the projector onto states symmetric under 
coordinate permutations and with parity ±1 under sign reversals, it is easy to show that 

A SC (L 2 (C)) = A+(L 2 (C)) ©A-(L 2 (C)) . (36) 

3.2. Basis of L 2 {C) 

Our next step is to construct suitable (non-orthogonal) base£]of A(L 2 (C) ® S) and 
A S c(i 2 (C)) on which H' <g> II and H', respectively, act as triangular matrices. The 
decompositions ([35 j) -([36 |) . and the fact that H' (8 1 and H' clearly commute with A ± 
and A^,, suggest that we first triangularize H' on L 2 (C). It should be noted that this 
problem is considerably harder than the corresponding one for the rational Dn model 



studied in Ref. 44 1, due to the fact that in the latter case C = = R w . However, in 
the present case C does not coincide with C' B ) = [— ^ , \ ] N , so that one cannot assume 
that the functions 

p(x) e 2inx , n 6 Z, N , (37) 



obtained from the basis of £ 2 ([— f , §] N ) found in Ref. [42j by setting b = b' = 0, are 
a basis of L 2 (C). In fact, it turns out that the set (|3"Tf is not complete in L 2 (C), and 
must therefore be supplemented by additional functions in order to obtain a basis. This 
peculiarity, which is absent in the rational case, lends an additional layer of complexity 
to the relation between the trigonometric Dn models and their BCn counterparts. 
In this subsection we shall prove that the functions 

pW(x) = p(x)e i £i< 2n ' +a > B ', n=(n Il ...,njf)6Z" 1 S € {0,1} , (38) 

form a Schauder basis of L 2 (C), leaving for the next subsection the proof that H' acts 
triangularly on the latter set when it is ordered appropriately. As we shall next discuss 
in more detail, the completeness of the functions ([38]) is essentially based on the fact that 
a complex exponential e lk x (k £ R) is periodic in C if and only if 

k= (2n 1 +6,...,2n N + S), rij G Z, 5g{0,1}. (39) 

Since C is not a hypercube, we need to define more precisely what it means for a 
function to be periodic in this set. To this end, let 

Fff = {x e C : x, + exj = e'tt} , 1 i < j < N , e, e' = ± , 

denote one of the 2N(N — 1) faces of C. If Ty is the translation along the vector 
e'7r(ei + ee^) perpendicular to the latter face (where {ei, . . . ,ej\r} is the canonical basis 



1 More precisely, a non-orthogonal basis of a (separable) Hilbert space f) is a Schauder basis of its 
underlying Banach space, i.e., a countable subset {vi : i £ N} C such that every element v G S) can 
be expressed in a unique way as J^Si c i v i> with £ C. In the rest of this section the term "basis" will 
often be used in this more general sense. 
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of R-^), then each T^j~ e clearly sends the corresponding face F^ to its opposite F^J~ e . 

Given a point x £ , we shall refer to T^~ e x as the point opposite to x in the 

face F^~ e . (Of course, a point lying on the intersection of k > 1 faces has k different 
opposites.) 

We shall say that a continuous function / : C — > C is periodic in C if 

/(x) = f{Ttp x) , Vxei^', l<i<j<7V, e,e' = ±. (40) 

In other words, / is periodic in C if it takes the same value on opposite points in any 
two faces of C. Since the coroots of the Dn root system with the normalization (fTTj) are 
the 2N(N - 1) vectors 

7r(±e l ± e,) , 1 sC i < j ^ N , 

the group X generated^ by the translations T" is the translation group corresponding 
to the Dm coroot lattice (i.e., the Z-linear span of the coroot vectors). As is well known, 
the semidirect product of 1 x 2B yields the affine Weyl group of Dn type 22J a [49j |. We 
shall define two points x, x' £ 'SL N to be equivalent, and shall write x ~ x', provided that 
x' = Tx for some T £ T. Note that Eq. (l40l) and the previous definition imply that a 
function / is periodic in C if and only if 

x,x' e dC, x ~ x' => /(x) = /(x') . 

We shall accordingly say that a continuous function / : M. N —> C is T-periodic if it 
satisfies 

/(x) = /(Tx), VT e X , Vx e R w . (41) 

Remark 3. Since T^t is a translation of 2ir in the direction of the vector ej, it 

follows that every T-periodic function is 27r-periodic in each coordinate. The converse, 
however, is not true in general. 

One of the main ingredients in the proof of the completeness of the set (|3"5)) in L 2 (C) 
is the fact that every continuous function / : C — > C periodic in C can be uniquely 
extended to a 1-periodic function / : — > C. This result is a direct consequence of the 
following fundamental facts: 

i) For each x £ M. N , there is a point x'eC such that x ~ x'. 

ii) Moreover, if x ~ x" £ C and x" ^ x', then both x' and x" lie on (at least) a face 
of C. 

The proof of these two statements is straightforward. Indeed, it is well known that 
A is a fundamental domain for the action of 2U a in R w , i.e., that for every x £ ~R N there 
is a unique & £ A and a suitable element R of 2B a such that x = i?a. Since 2U a is the 
semidirect product of its subgroups % and 2U, we can write R = TW, with T £ % and 



Ty' for i < k < j, the group X is generated just by the translations 
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2 In fact, since T+' T T.+' 

rp+,± 

ij 



W £ 20. This shows that x is equivalent to x' = Wa. 6 W(A). Since the elements of 20 
are homeomorphisms, we have W(A) — 20(A) = C, which proves the first statement. 

As to the second one, suppose next that x is equivalent to two different points x' and 
x" of C. It follows that x' ~ x" or, equivalently, 

x 'l = x 'i + kiir , ki £ Z , 1 ^ i ^ N ; k\ + ■ ■ ■ + kpj £ 2Z . 

Since both x' and x" belong to C C [— 7r,7r] , the integers fcj can only take the values 
0, ±1,±2. Suppose, first, that one of these integers is equal to ±2. Without loss of 
generality, we may assume that k\ — 2. Since x'i,x" £ [— 7r, 7r], this is only possible 
if x[ = —it = —x'l, from which it follows (taking into account that x',x" £ C) that 
x\ = x'l = for i > 1. Thus in this case x' belongs to the 2(N — 1) faces F^'~ (with 
2 < i < N), and x" to the opposite faces F^' + . (Note, however, that in this case x" is 
not the opposite point of x' in any of these 2(N — 1) faces.) 

Assume now that |fcj| < 2 for all i. Since not all the integers fcj can be zero by 
hypothesis, and k\ + • • • + k^ must be even, there are two of these integers, which 
w.l.o.g. we can take as k\ and fe, such that k\ — ek^ = e' (with e,e' = ±). From the 
equalities x'{ + ex'2 = x[ + ex' 2 + 2e'ir and the fact that x[ + ex' 2 , x'{ + ex' 2 ' £ [— 7r, n], it 
follows that x'{ + ex'2 = ~( x 'i + ex 2) — e,7r i an d therefore x' and x" lie on (at least) the 
opposite faces and F^ 2 ~ e of C, respectively. This completes the proof of the second 
statement. 

Remark 4. The previous result implies that iV-dimensional Euclidean space can be 
tiled with translations of the set C along the Dm (co)root lattice. For N — 3, this is the 
well-known tessellation of R 3 with rhombic dodecahedra [H(| • 

We are now ready to prove that the set (|3"5)) is a (non-orthogonal) basis of L 2 (C). In 
fact, it suffices to show that the exponentials 

e i ^(2"i+*)^ j n=( ni ,...,n N )€Z, S £ {0, 1} (42) 

are themselves a basis of L 2 (C). Indeed, if this is the case then any complex- valued 
function / £ Cq(C) continuous in C and with compact support in C can be represented 
in the form 

/(x)= ]T c^e^ 2 "^^, (43) 

neZ N , 5£{0,1} 

where the coefficients c n j € C are by hypothesis uniquely determined by /. Our claim 
follows immediately from the fact that the function f / p is also in Co(C), and that the 
latter set is of course dense in L 2 (C). 

The fact the exponentials ((42]) form a basis of L 2 (C) is essentially a consequence of 
the fact that the momenta (|39p are the elements of the dual (or "reciprocal" , in a more 
physical terminology) lattice of the Dm coroot lattice. However, for completeness' sake 
we shall next provide an elementary proof of this fact. We first note that, since the set 
P(C) of complex- valued continuous functions periodic in C contains the dense set Co(C), 
to prove that (j42|) is a basis of L 2 (C) we need only show that every / £ P(C) can be 
uniquely represented by a Fourier series of the form (|43)l . Let, then, / : C — » C be a 
continuous function periodic in C, and denote by / : — > C its T-periodic extension. 
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Since the function / is 2-7r-periodic in each coordinate (cf. Remark [3]), it can be developed 
in terms of the L 2 ([-ir, ir] N ) Fourier basis e ikx (k e Z N ) as 

/»= £ a k c ikx , (44) 

keZ« 

with a k € C uniquely determined by /. Imposing that / satisfy Eq. (|41j) when T is one 
of the generators T?j of T we easily obtain 

a k (l-e i7r£ ' (fe!+£fe ^) =0, VkeZ^, 1 i < j TV , e,e' = ±. 

Hence a k vanishes unless 

fci ± fcj ee kfj G 2Z , 1 sC i < j sC N , 

i.e., unless all the integers fej have the same parity. Setting fc, = 2rii + S (with (5 € {0, 1}), 
&k = Cn.5, and substituting into Eq. f)44[) we obtain 

nez N , <5e{o,i} 

where the equality should be understood in the sense of L 2 ([—tt,tt] n ). Restricting to 
L 2 (C) (which is allowed, since C C [— 7r,7r] ), and recalling that all the coefficients a k 
are uniquely determined by /, we obtain the desired result. 

Remark 5. In fact, since the translations of C along the Dm (co)root lattice are a tiling 
of R N , the results of [5l| imply that the functions (|4"2l are mutually orthogonal. Of 
course, this does not imply that the functions ([55)1 are themselves orthogonal, due to the 
presence of the factor p(x). 

3.3. Triangularization of H' 

We shall next endow the set (|38p with a suitable order such that the action of H' on 
the resulting basis is triangular. Note, first of all, that 

L 2 (C) = SjW®& 1 \ (45) 

where is the closure of the subspace spanned by the basis functions ip^ with n G Z w . 
We will show that H' leaves invariant each of the subspaces JftW, so that we need only 
order each subbasis {ip^ } neZ « in such a way that H' is represented by a triangular 
matrix in fj( 6 '. To this end, given a multiindex p ee (p 1; . . , , J»jv) € % N we define 

[p] = (biil>---,bwl) , with \Ph\ > ■•• > bijvl- 

If p' e Z w is another multiindex, we shall write p -< p' provided that the first non- 
vanishing component of [p] — [p'] is negative. The basis functions {Vn } neZ jv should 
then be ordered in any way such that tp^ precedes (p^) whenever u -< u', where 

u ee {2m +6,...,2n N + S), (46) 
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and similarly for u' . For instance, wjj j \ must precede ip 



(%_ 3 ^(31 0) s h ou ld follow 



¥?(2 _3_i)> while the relative precedence of -1) an d V(i 3 -2) can ^ e arbitrarily 

assigned. 

In order to compute the action of H' on the basis functions (|38[) , it is convenient to 
introduce the Dunkl operators of Dn type 

Jk =id Xk +aJ2 [0--icotx^)K k i + (l-icotx^)K k i -2a^X fc; , (47) 



l^k 



Kk 



with k — 1, . . . ,N, obtained from their BCn counterparts [42| by setting b = b' = 0. 
Note that the set C is invariant under all the generators Kij, Ki of the Weyl group 
of BCn type 2U( B \ and hence the operator J k is actually defined in a suitable dense 
subspace of L 2 (C). It can be shown that 



(48) 



so that the action of H 1 on the basis (|38|) can be easily inferred from that of the Dunkl 
operators (|47|) . In the following discussion, we shall label the basis functions ipn simply 
by <p v , with v given by (|4"fJ|) . As in Ref. [521 ] . we shall begin by considering the action of 
J/j on a basis functions ip„ with 1/ nonnegative and nonincreasing. For such a multiindex, 
we shall use the notation 



#(s) = card{i : vi = s} . 

with £(s) = +00 if vi 7^ s for all i = 1, . . 
#(2) = 3 and £(2) = 4. 

We shall next prove the formula 



£(s) = min{« : = s} . 
, N. For instance, if v = ( 



,6,6,2,2,2) then 



Jk^v = At,,fc ip v + 



E 



C„,fc <Pv> 



v' -ve(2Z) N , v'-< v 



(49) 



where k e 



and 



-1* + 2a(2£(u k ) + #{v k ) -k-N-1), v k > 
2a(iV - fc) , j/fc = , 



(50) 



which will play a fundamental role in the sequel. To begin with, a lengthy but straight- 
forward calculation yields 



J ^ = -v k -2a(N-l) + 2aJ2 — 



1 



+ 2 «E 

i>fe 
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2«E 

3+h 



n2—Vj—u h 
Pjk 



/3 2 



j A- 



(51) 



where 

Oij k — Zj Zfc j fijk — ZjZk ) Zj = G ^ . 

Consider now the first sum in Eq. (1511) . Since j < k, by hypothesis Vj ^ z>%. If Vj = Vk, 
the j-th term in this sum clearly vanishes. On the other hand, if Vj > Vk we have 

<Pu % —, = <P» + E z7 2r 4 r V», (52) 

where the last sum only appears if Vj — v k > 2. In this case the multiindices u' of the 
monomials in the summation symbol in Eq. (|52j) are of the form 



(vi, . . . - 2r, . . . , v k + 2r, . . . , v N ) , r = 1, . . . , ~(vj -v k )-l. 



and hence v' — v € (2Z) . Moreover, we have < maxjiA,- — 2r, v k + 2r} < i/j for all 
r = 1, . . . , — j/fc) — 1, so that i/ -< i/. Thus, the first sum in ([51]) contributes to A„.fc 
the quantity 

2acard{j < k : Vj > v k } = 2a(l(v k ) - l) . (53) 

It may be likewise verified that the multiindices z/ corresponding to the monomials 
arising from the second sum in Eq. (1511) cither coincide with v or satisfy v' -< i/, and 
that this sum yields the following contribution to \ v y. 

2acard{j >k:v J - v k } = 2a{t(v k ) + #(i/ fc ) - fc - l) . (54) 

Consider, finally, the j-th term of the last sum in Eq. (|5ip . This term is equal to 1 when 
Vj = v k = 0, while for Vj + v k ^ 2 we have 



Hence when ja,- + z/fc ^ 2 the multiindices corresponding to the basis functions arising 
from the last sum in Eq. (I51[) are of the form 



v' = (z/i, . . . , Vj - 2r, . . . , v k - 2r, . . . , v N ) , r = 1, . . . , -(vj + v k ) - 1 , 

and thus v' — u g (2%) N . Furthermore, since 

-i/fc + 2 ^ < I/,- - 2 , -i/j ■ + 2 < i£ < z/ fe - 2 , 

we have max{|z/|, \v' k \} < maxjzA,-, v k }, so that again v' -< v. The contribution of the 
last sum in Eq. ((511) to is therefore equal to 

2a(#(0)- 1)^,0. (55) 

Adding Eqs. (I5"3l - (j55|) to the first two terms in the RHS of Eq. (|5"Tj) . and taking into 
account that 1(0) + #(0) = N + 1, we easily obtain Eq. (|5D|) for A^fe. 
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Equation ([49]) does not hold in general if v does not belong to [Z w ] , so that Eq. ((50|) 
does not yield the spectrum of the Dunkl operators J k . On the other hand, for the 
purposes of computing the spectrum of H' we shall only need the following weaker result: 
if v G 1i N is a multiindex all of whose components have the same parity, then 

Jk<Pv = l„,k <Pw (56) 

u'ez N 

v'-ve(2Z) N , [v'\^[v\ 

for some complex constants 7^ k . Indeed, if v is as above, there is an element W belonging 
to the Weyl group of BCn type such that ip v — Wip[ u ]. Proceeding as in Ref. [H^ . 

it is easy to show that 

2 N iV! 

[J k ,W] = C 3* W 3> c jk eR, 
where W (B *> = {Wj : j = 1, . . . , 2^!}. We thus have 

2 N N\ 

J k ip v = W(Jk(P[u]) + X! C Jk W 3 

and Eq. follows immediately from (l4"9l and the fact that the partial ordering -< and 
the parity of the components are invariant under the action of 2Jj( B ) . 

From the previous results it is relatively straightforward to compute the spectrum of 
H' . More precisely, we shall next show that the action of H' on each Schauder subbasis 
{^n } n g Z » , ordered as explained above, is upper triangular: 

H'<pV> = E<?<p<$ + ]T , v k ee 2n k + S , v' k ee 2n' k + 5 , (57) 



where c£/ n € C and 



E^=^2([^] k + 2a(N-k)) 2 . (58) 



Indeed, suppose first that the multiindex v in Eq. (|57p is nonncgativc and nonincreasing. 
Applying Jk to both sides of Eq. (|4"9"|) we obtain 

Jk<Pu = >i,k<Pv + X »' k C ",k <Pl>> + Y C ";k J W"' ■ 

By Eq. e d56l the last sum is a linear combination of basis functions ip„n with v" -< u 
and v" — v G (2Z) N . Therefore we can write 

J k cp v = Xl jk (f v + Y b »,k <Pv' > 
i/'-i/e(2Z) N 

V -< V 
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with k E C. Summing over k and using the identity (|48|) we obtain 



#v = (E A ^K + E (E<*V- ( 59 ) 



i/-^e(2Z) f 



Suppose, next, that ^ [Z]^, and let W G 2H^ B ^ be such that y>„ = Wipy v y Since 
ff' is obtained from its BCn counterpart in Ref. 42] by setting b = b' — 0, and the 
latter operator commutes with all the elements of W^ B \ it follows that [H ! , W] = 0. By 
Eq. (|59l) applied to we have 



ffV^=w-^v M = (E A kfe)^+ E (E 6 m,*K^- 

which establishes ((57)) with 



i/'-He(2Z)' 



All that remains to be proven is Eq. (|55|) for the eigenvalue i?i' 5 ' ) . To this end, let p = [u] 
and suppose that pk-i > Pk = ■ • • = Pfe+r > Pfe+r+i 0, so that £(j>k+j) — k and 
#(.Pfc+j) = r + 1 for j = 0, ... , r. Since 

Ap,fe+j = -Pfc+j + 2a(fc + r-j - N) = -p k + r -j + 2a(fc + r - j - N), j = 0, . . . ,r , 

we have 

E A L = Efe + 2a ( iV -^) 2 - ( 61 ) 

j=fc i=fc 

If, on the other hand, pk-i > Pk — • • ■ = Pn — 0, the analog of Eq. (|rJT|) follows directly 
from (f5U|) . This completes the proof of Eq. 



5.^. Triangularization of H and H sc 

Using the results of the previous subsection, it is a straightforward matter to trian- 
gularize H and -ff sc . Indeed, by the results in Subsection 13. 1[ this problem is equivalent 
to the triangularization of the extensions H and H sc acting on their respective Hilbert 
spaces ft = A(L 2 (C) ® S) and Sj sc = A SC (L 2 (C)) 

Let us start with the operator H. By Eqs. (|3"5j) and P5|) . its Hilbert space can be 
decomposed as the direct sum 

^=0 A e (i^®£). (62) 

<5=0,1 

Let /(x) be in the domain of H', and let |s) € E denote an arbitrary spin state. Since 
if coincides with H 1 ® 1 on 55, and the latter operator commutes with A (indeed, it 
commutes with all the elements of W^ B \ and hence of 20), we have 



ff[A e (/(x)| s ))] =A*[(Jf'/(x))| S >] 
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(63) 



As we saw in the previous subsection, H' preserves the subspaces Sj^ , which by the latter 
equation implies that each of the four subspaces A e ^^ (g> S) is invariant under H. We 
shall next verify that H acts triangularly on a (non-orthogonal) basis of A e ($j( s ' (g) S) of 
the form 

^ s (x)=A e (^)(x)|s)), (64) 

ordered in such a way that ip^ e s precedes tp^fa' whenever v -< v' (with 1/ defined in (|4"6")l . 
and similarly i/). The spin functions (|64p are obviously a complete set (since the func- 
tions |38|) are a basis of L 2 (C)), but their linear independence is not assured unless 
we impose suitable restrictions on the quantum numbers n and s. More precisely, the 
states l|64p are a (non-orthogonal) basis of the Hilbert space A e ($j( s ' ® S) provided that 
the quantum numbers n 6 Z w and s satisfy the following conditions: 

i) n\ > • • • ^ n^v 5^ (65a) 

ii) s 2 ; > Sj whenever m = Uj and i < j (65b) 

iii) If S = rii = then Si ^ for e = 1, while Sj > for e = — 1. (65c) 
Indeed, since 

A e (KiSi) = eA e , A'(KijS i:i ) = -A e , 

acting with suitable operators XiSi and KijSij on a spin function <£>n (x)|s) with arbi- 
trary n G Z w and s one can easily show that the corresponding state i/i* ,e s is either zero 
or proportional to a state satisfying the above conditions. (Note, in this respect, 
that a state (IMl) with 5 = = = is symmetric under (xi,Si) — > (— Xj,— Sj), and 
must therefore vanish identically if e = —1.) This shows that the states with n £ 1, N 
and s satisfying the above conditions are complete. Their linear independence is easily 
checked. 

Remark 6. Conditions i) — ii) above are identical to the corresponding ones for the spin 



Calogero model of Dn type studied in Ref. [44[. As to the third one, the key difference 



is that in the present case the action of a coordinate sign reversing operator Ki on a 
state <£n^(x)|s) no longer produces a state with the same quantum number n (up to a 
constant factor) unless 5 — n% = 0. 

Remark 7. Since the functions <^ 0) with n <= Z N form a basis of L 2 {C {B) ) C L 2 (C) 
(cf. Ref. [42|]), it follows that each subspace A e (Sy^ g3 S) properly contains the Hilbert 
space A e (L 2 (C( B )) ® S) of the Sutherland spin model of BCn type ([T]) with chirality e. 
Note, however, that the other sector A + (ijWcgiE) © A~ f.£jW(g)£) of has no counterpart 
in the BCn model. Thus, in contrast with the rational case [44|, the Hilbert space of the 
Dn model is larger than the direct sum of the Hilbert spaces of its BCn counterparts 
with both chiralities. 

Let us now examine the action of the operator H on the basis of A e (Sj^ ® S) given 
by Eqs. (|64[) - (f^5|) . It is easy to show that 

= < £ s V4',s + E c »v, ns . (66) 

n' ,s' 
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where c^, s , ns £ C and 

E 5 ^ s = J2( 2n k + S + 2a(N -k)Y . (67) 

k 



Indeed, from Eqs. (j57l) - (j58|) and the identity (|63| one immediately obtains 

= + ]T <#> n . (68) 

Although the quantum numbers (n',s) appearing in the RHS of this equation do not 
necessarily satisfy conditions (p5|) , there is an element W £ such that (Wri, Ws) = 

(n",s") do satisfy these conditions. Since the corresponding state ip^f, s „ differs from 

tp^s ^ a ^ mos t an overall sign, and \v"] 
that we can rewrite (|6"8")l in the form (|6"6")) . 



V'n^s ^ a t mos t an overall sign, and [v"\ = [u'\ -< [v] implies that v" < v. it is clear 



From Eq. (166[) it follows that the operator H acts triangularly on the (non-orthogonal) 
basis of A e (Sj( 5 ' ® E) in Eqs. (jMl) - (l6"5l) . ordered in such a way that ip^\ precedes ip^f B i 
whenever u -< i/. Moreover, the eigenvalues of the restriction of H to A £ (f)( 5 ) <g) E) are 
given by Eq. (|67[) . with n S Z w and s satisfying conditions (p5|) . 



Remark 8. Since the numerical value of the eigenvalue (I57|) does not depend on s or 
e, for any multiindex n £ [Z w ] the corresponding eigenvalue E^ e s has an associated 
degeneracy 

d 5 n = tg+ + dj- , (69) 

where d* ,€ is the number of basic spin states |s) satisfying conditions (|65p for given e and 
5. These spin degeneracy factors will be computed below, when we discuss the partition 
function of this model. 

The spectrum of the scalar Hamiltonian H sc can be computed in a similar way by 
exploiting the fact that it coincides with H' in its Hilbert space fj sc , which by Eqs. |36|) 
and is given by 

£ sc = Al c {SjW). (70) 



8=0.1 



Due to the identity 



# sc (A| c /(x))=A s e c (ff'/(x)) 



it is immediate to show that each of the four subspaces A^i}^) is invariant under H sc . 
Just as in Section EOT) . it can be verified that the functions 



^ e (x)=A| c (^)(x)), (71) 

where n £ 7L N and 

ni>--->n^>i(l-e)(l-5), (72) 

are a Schauder basis of k\ c {%^^. (The last inequality is due to the fact that if 6 = tin = 
the function Vn'~ is symmetric under xn — > —xpj, and therefore vanishes identically). 
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Proceeding as above, it is straightforward to show that if we order the basis (|71[) - ((72|) 
so that i/A e precedes tp^ whenever u -< u', the operator H sc acts triangularly on it, 
with eigenvalues given by the RHS of Eq. (p7|) . Of course, due to the absence of 
internal degrees of freedom, in this case the degeneracy factors d^ e are equal to one for 
all quantum numbers n, e = ±1, and 5 = 0, 1. 

Remark 9. It is well-known [22j that the eigenfunctions of the scalar BCn Sutherland 
model are of the form 

p(x) H| sin Oleosa;/ ■ P k {y) , (73) 

i 

where -Pfc(y) is a symmetric polynomial in the variables yi = sin 2 Xi (i = 1, . . . , N). The 
polynomials P k , which can be regarded as multivariate generalizations of the classical 
Jacobi polynomials, are orthogonal in the hypercube [0, 1]^ with respect to the weight 
function 

^ (B) (y) = n \i* - ytf a ■ IT y"~^ 1 - yJ b '~ h ( 74 ) 

i<j i 

(cf. Eq. (2.17) of Ref. (22[). In our case, from the identities 

A+ e 1 ^ = A sc J] cos(v k x k ) , A- e 1 ^ v " x " = i N A sc ]J sin(v k x k ) , 

k k 

it is straightforward to show that the (orthonormalized) eigenfunctions of H sc in each of 
the invariant subspaces Ag C (^^^) are of the form 

P(x) J| | sin 2xi\^\ cos n \ 5e ■ P*' e (y) , (75) 

i 

where P k e {y) is a polynomial in the variables yi = sin 2 Xj, symmetric under permutations. 
From the discussion in Section 13.11 it follows that the restrictions of the functions (l75l) 



to the open set A are a complete set of eigenfunctions of the scalar Sutherland model of 
Dn type H sc . They are also orthogonal in the latter set, on account of their symmetry 
under coordinate permutations and sign changes. This is easily seen to imply that the 
polynomials Pfj e are orthogonal in the hypercube [0, 1]^ with respect to the weight 

^■ £ (y) = I]> - yj\ 2a • n^ r *(i - VrY (S - h) ■ (76) 

i<j i 

In view of Eqs. (j73l) - ([71|) and (|75p - ([75|) . it is clear that the three orthogonal polyno- 
mial families {P k ' e : fc € N} with (5, e) = (0, — 1), (1, ±1) are not limiting cases of the 



multivariate Jacobi polynomials studied by Baker and Forrester 22]. The analysis of 
the properties of these new orthogonal polynomials, and their relations with their BCn 
counterparts, could lead to interesting new developments in the field of multivariate 
orthogonal polynomials. 
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4. Partition function of the spin chain 



The purpose of this section is to evaluate in closed form the partition function of the 
Haldane-Shastry spin chain of Dm type (fTB"]). Our starting point is the freezing trick 
relation (|17p . which can be equivalcntly written as 

£ = ii m ^ ^ . (77) 

a^oc 4a 

This formula expresses each eigenvalue £j of the chain (JT5]) in terms of a certain eigenvalue 
Eij of the spin Sutherland model of Dn type (O and a corresponding eigenvalue Ef c of the 
scalar model (TT6]) . In practice, the fact that the eigenvalues J3y and Ef c are obviously not 
independent makes it impossible to use Eq. (1771) to completely determine the spectrum 
of the chain (|T3"f in terms of the spectra of the Hamiltonians H and H sc computed in the 
previous section (cf. Eq. (|67|) V The key idea behind the freezing trick method introduced 
by Polychronakos [36| is to use Eq. ((77)) , or rather the equivalent relation fTT)) , to directly 
compute the chain's partition function. Indeed, the latter equation immediately yields 
the identity 

Z (T) = lim Z{ f aT \ , (78) 

expressing the chain's partition function Z in terms of the partition functions Z and Z sc 
of the Hamiltonians H and H sc . 

Remark 10. Equation (|77[) can be used to obtain nontrivial qualitative information on 
the spectrum of the chain (|13[) . For instance, from the fact that the numerical values 
of the energies of both Hamiltonians H and H sc are given by the RHS of Eq. (p7|) and 
Eq. (|77p it easily follows that all the energies of the spin chain ([T3"| are integers. 

In the rest of this section, we shall compute the large a limits of Z(4aT) and Z sc (4aT) 
using Eq. (|67[) for the spectrum of H and H sc , thereby obtaining an exact expression for 
Z via Eq. (|78l) . Before doing so, it is convenient to subtract from the spectra of H and 
H sc the constant term 

E = 4a 2 ^(7V - k) 2 = | a 2 iV(7V - 1)(2N - 1) , 
k 3 

which is of course irrelevant for the purposes of computing Z. The rationale behind this 
normalization is the fact that, by Eq. (|57|) . the eigenvalues of H and _ff sc become 0(a) 
for a — > oo, so that the limits of Z(AaT) and Z sc (4aT) exist separately. 

Let us start with the partition function of Hamiltonian H of the D^-type spin Suther- 
land model (JSJ). With the normalization of the energies explained above, the spectrum 
of this model satisfies 

£n's = 4a^(2n fe + 5)(N - k) + 0(1) , (79) 

k 

and hence its partition function is given by 

d s n e q> , g-e- 1 /^). (80) 

ni -5? ■ ■ ■ ^2 n /v 5?0 
e=±, 5=0,1 
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As mentioned in Remark the degeneracy factor c?* ,£ is equal to the number of spin 
states |s) satisfying conditions (j65|) for given e = ±1 and <5 = 0, 1. Writing the quantum 
number n in the form 



n= (p 1 ,...,px,...,p r ,...,p r ), 
and using conditions (|65b|) and (|65cp we have 



n 



n 



with 



Let us now define 



1 / 
= - 2 {m 

Z< 4 >(T) 



r(m)) , 



Pi > • • • > p r ^ 0, 

5 = p r = ; 
otherwise, 



7r(m) = m (mod 2) . 

E(2ni+<5)(W-i) 



(81) 



(82) 



(83) 



e=rt 



so that 



lim Z(4aT) = Z<°>(T) + Z {1 \T) . 



(84) 



The function Z^> (T) can be easily expressed in terms of the partition functions Z^ of 
two spin Sutherland models of BCn type ([1]) with opposite chiralities e = ±. Indeed, 
when (5 = Eq. (|82p coincides with Eq. (51) in Ref. 42[ for the degeneracy factor of 
the -BCV-type spin Sutherland model with chirality e. Likewise, Eq. (|6T|) with <5 = 
is obtained from the analogous formula for the energies of the BCn Hamiltonian ([1]) in 
Eq. (24) of the latter reference by setting /3 = /3' = 0, and the same is true for Eq. (1751) . 
We thus have 

Z {Q) (T) = lim (zi B) (4oT) + zi B) (4aT)) . (85) 

a-s-oc V "f / p—pi—Q 

Using Eq. (53) from Ref. [42| we obtain the explicit expression 



z<°>(t)= ^ 

(fei,...,fc r )ePjv 



m_ 



K r 



n 



-,A\ 



1 - <z A '* 



, (86) 



where Vn is the set of partitions of the integer N (taking order into account), and 



K i = k i (2N -l-h) 
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j=i 



(87) 



Note that, since k\ + ■ ■ ■ + k r = N, the integers fc, are in the range 1, . . . , N. 
On the other hand, from Eq. (|82p it easily follows that 



Proceeding as in Ref. [42| we easily obtain 



zW(d = 2^-) 2 a - **r i n ( D • n 

(fci,...,fc r )ePjv <=i v 4/ »=i 

The partition function Z sc (4aT) of the scalar Hamiltonian ([TBI) is also easily evaluated 
in the limit a — > oo, since in this limit its spectrum (with the normalization discussed 
above) is still given by the RHS of Eq. (|79[) . Using Eq. (|72l) . and taking into account 
that in this case <i* ,£ = 1, we have 

E(2n*+l)(JV-i) E2n,(7V-i) 

lim z sc (4aT) = 2 > g i + > 

a— >oo * — » * — » 

£2fn(JV-i) 

+ E 

= (2g^(JV-i) + 1) ^ r + £ q i 

7ii ^ ■ ■ ■ ^ n iv ^ 

The last sum is easily recognized as the a — > oo limit of the partition function Z\ c ' (4aT) 
of the scalar Sutherland Hamiltonian of BCn type 

= 



with /3 = j8' = 0. Using Eq. (49) in Ref. [42j we thus obtain 
lim Z sc (4aT) = (l + q^^- 1 )) 2 lim Zi c B )(4aT)| 

a— >-oo a^-oo 'P— P — u 

= (1 + g w- 1 )) 2 rj (i - g^tf-i-f))- 1 . (89) 

The partition function of the Haldane-Shastry spin chain of -Djv type (|13[) is easily 
computed by inserting Eqs. (J84J), (|86|>. (|88)l and flSS) into the freezing trick identity _d78l) . 
In order to simplify the resulting expression, we define N — r integers k[ < ■ ■ ■ < k' N _ r 
in the range 1, . . . , N — 1 by 

{k[,...,k' N _ r } = {l,...,N-l}-{h,---,k r -i}, 

and set 

Ki = k' i (2N-l-k' i ). (90) 
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Using this notation, the partition function of the chain (1131) can be written as 

„ i r-l 



Z(T)=(l + q i^-V)- 2 E II u I 

(k 1 ,...,k r )£V N i=l 



1-9 



A'r 



JV-r 



(91) 



Taking into account that k r — ki — N, so that K r = N(N — 1) by Eq. (|57|) . we finally 

i=l 

obtain the more compact expression 



(k lt ...,k r )eV N i=l 



TO- 



2| ™] g *"("-D 



AT— r 



-A" 



(92) 



Remark 11. From 0, Eq. (53)] and Eq. (ItJTI) we easily obtain the identity 



Z(T) = (1 + ^ W ( w -D)- 2 [fzf >(T) + Z^(T)) + 2q? N ( N -VQ N (T) 



where 



Qim= e n 



(fei 



,k r )eV, i=l 



r - i l-r 
i=l 



8^0 



(93) 



and the integers Ki, K[ are defined by Eqs. ([57)1 and (|§U1) for all Z. Thus, unlike what 
happens in the rational case (cf. [3, Eq. (46)]), it does not seem possible to express in a 
simple way the partition function Z(T) exclusively in terms of its BCn counterparts Z ± . 
Note also that the function Qat(T) has the same structure as the partition function of the 
ordinary (Ajy-type) Haldane-Shastry chain 37j . the only difference being the "dispersion 
relation" defining the quantities Ki and K[ in terms of ki and k[. 

As mentioned in Remark 1101 the eigenvalues of the spin chain (1131) are integers, and 
they are nonnegative on account of the nonnegative character of the operators 1 + Sy 
and 1 + Sij. Thus, the partition function Z(T) should be a polynomial in q, a fact which 
is not apparent from Eq. (|92p . In order to ascertain this fact, consider first a partition 
(k\, . . . , k r ) 6 Vn with k r = 1. In this case the term in curly brackets in Eq. (|92p reduces 
to m(l + g Ar ( Ar - 1 )/ 2 ), and k r -i = N-k r =N-l implies that 



{k[,...,k' N _ r }={l,...,N-2}-{k 1 ,.- 



Hence the contribution to Z(T) of the partitions with k r = 1 is given by to q N ( N ~ 1 " > Qn-i(T). 
Consider next a partition (k±,...,k r ) £ Vn such that k r = I > 1. In this case 
(fci, . . . , fc r _i) is a partition of N — I, and /c r _i — N — I implies that 



J - 1 



(94) 



2(i 



and hence 



K' N ^_ r+1 = (N-j)(N + j-l) 



j = l,. ..,1-1, 



Note, in particular, that K' N — N(N — 1), so that 

(1 + gW- 1 )) (1 - = i _ q iN(N-l) 

Taking into account that, by Eq. ((94| . 

{fci,...,fc^_ r+1 } = {l,...,N-l-l} - {h,--.,k r -2}, 

it is immediate to verify that the contribution to 2(T) of the partitions with k r = I ^ 2 
is given by 



(1 



1-2 

g iiY(JV-i)yjv-()(jv+;-i) j| ^ 



2| ■|^ g ^(AT-l) + 



(JV-i-l)(AT+i)> 



A i 



ra_ 



1 



iA^(^'-l)^ 



2»-;(T), 



with Qo = 1. Thus the partition function (|92l) can be expressed af0 
m q 



Z(T) = mq N ^-^Q N -i(T) 



1 



±JV(JV-1)"\ 



min(m, TV) 



E 



(jv-o(jv+j-i) 



nc 1 



2=2 



V l 



I 



1-9 



9 

iiV(JV-l) 



) Qat-^T), (95) 



where the RHS is clearly a polynomial in q on account of Eq. (|93[) . This remarkable 
formula is one of the main results in the paper. 

Remark 12. From Eqs. (l!?2"j) or (p?5)) it is apparent that the partition function of the 
Dn chain (|13[) has a much more complex structure than its BCn counterpart with 
P = j3' = 0, cf. 42J, Eq. (53)]. In particular, while for the BCn chain one can find [53| a 
description of the spectrum in terms of a suitable generalization of Haldane's motifs [541 ] . 
it is not clear how to implement such a description for the present chain. Note that, for 
HS chains of An type, the existence of such a description is the key ingredient in the 
proof of the Gaussian character of their level density [55] when the number of sites tends 
to infinity, which is of importance in the context of quantum chaos. 



5. Concluding remarks 

As mentioned in the Introduction, reductions of the BCn Calogero and Sutherland 
models obtained by setting suitable coupling constants to zero have been largely ignored 



3 Note that the terms with I > m in the previous equation vanish identically due to the binomial 
coefficients. This is in fact a consequence of conditions l !65bl l and Il65cl) . cf. Eq. 1821 . 
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in the extensive literature devoted to these models. This is probably due to the fact that 
these reductions were mostly regarded as trivial limits of the above models. In a previous 
paper [44J , we showed that this is not case by studying the Dn reduction of the (spin) 
BCn Calogero model. Moreover, the spin chain of Haldane-Shastry type associated with 
this reduction was also seen to differ from its BCn counterpart even more markedly, 
essentially due to the nontrivial nature of Polychronakos's "freezing trick" . The aim of 
the present paper is to perform a comprehensive study of the Dn reduction of the BCn 
spin Sutherland model and its associated spin chain. 

A significant part of the paper is devoted to the exact computation of the spectrum 
of the dynamical spin model © and its scalar version (|16l) . We have first provided a 
rigorous proof of the equivalence of these models to their extended versions H and H sc 
defined on the Weyl-invariant configuration space C. The latter set, which turns out to be 
the iV-dimensional generalization (1271) of a rhombic dodecahedron, is more complicated 
in nature than its BCn counterpart (a hypercube). The motivation for constructing the 
extended operators H and H sc is the fact that on their natural domains they essentially 
coincide with the restriction of a simpler auxiliary operator iJ', which can be expressed 
as a sum of squares of a suitable set of Dunkl operators of Dn type. In particular, 
from the spectrum of H' it is not difficult to deduce those of H and H sc , and hence 
of H and H sc . In order to compute the spectrum of H', we have next constructed an 
appropriate (non-orthogonal) basis of the Hilbert space L 2 (C) where this operator acts. 
This is indeed the key difference with the rational case, for which this step is trivial due 
to the fact that the (extended) configuratio n sp aces of both the Dn and BCn models is 
M. N . Using a method similar to that of Ref. [52j, we have shown that H' acts triangularly 
on the above basis when ordered appropriately. Finally, we have shown how to construct 
from the latter basis a (non-orthogonal) basis of the Hilbert spaces of H and H sc where 
their action is also upper triangular. In this way we have computed in closed form the 
spectra of the spin Sutherland model of Dn type (JSJ) and its scalar version (TlTJl) . 

The second main result in the paper is the exact computation of the partition function 
of the Haldane-Shastry spin chain of Dn type (|13|) obtained from the spin dynamical 
model ([5]) by means of Polychronakos's freezing trick. The latter chain, as is appar- 
ent from Eq. (|25p . cannot be obtained from its BCn counterpart by taking the limit 
(/?,/?') — > due to the presence of an "impurity" term at both endpoints. Our starting 
point is the fundamental relation ([75)1 . expressing the chain's partition function as the 
large coupling constant limit of the quotient between the partition functions of the corre- 
sponding spin dynamical model H and its scalar version H sc . Using the above mentioned 
results for the spectra of these models, we have been able to evaluate this limit, thereby 
obtaining Eq. (I92[) for the chain's partition function. In contrast with the rational case 
(cf. Remark I lip , the partition function is not expressed in a simple way in terms of 
its BCn counterparts, since it also involves the partition function of the original (type 
A) HS chain with a slightly different dispersion relation. We have further simplified 
Eq. (|92[) for the partition function, showing how to write it explicitly as a polynomial 
in q = e _1 /( fcBT ) (see Eq. (j9"5)) '). as should be the case for a finite system. In fact, this 
simplified formula turns out to be quite efficient for the numerical computation of the 
chain's spectrum, making it possible to perform a statistical analysis of the spectrum 
when the number of particles becomes very large. It would be worthwhile to carry out 
such a stud y, a nd com par e its results with the corresponding ones for other spin chains 
of HS type [37l [H 0, . 
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The results of this paper suggest a number of further developments that we shall now 
discuss. In the first place, we have shown that the Djy reduction of the standard Suther- 
land model of BCn type gives rise to an interesting new solvable model that had been 
previously overlooked. In fact, there are several additional reductions that could be con- 
sidered, like e.g. those associated with the Bn and Cn root systems, or even more general 
ones, like the 6 = reduction of the Sutherland model ((TJ) . It could also be of interest 
to consider similar reductions of the comparatively less studied hyperbolic Sutherland 



model of BCn type [52j. As we have also shown in this work, these reductions can 



potentially yield unexpected results in other fields, as for instance the remarkable tiling 
of M. N with the TV-dimensional generalization of the rhombic dodecahedron uncovered in 
Section [3] (see Remark 0]). 

The work presented in this paper has direct implications in the field of multivariate 
orthogonal polynomials. More precisely (see Remark [S]), the eigenfunctions of the Dn 
reduction of the scalar Sutherland model yield new families of multivariate orthogonal 
polynomials that cannot be obtained as straightforward limits of the generalized Jacobi 
polynomials associated with the BCn Sutherland model [22j ■ It is to be expected that the 
additional reductions mentioned above could lead to similar new families of orthogonal 
polynomials. 

An important aspect of spin chains of Haldane-Shastry type that has not been dealt 
with in this paper is their intcgrability, which for the original HS chain of An type was 
established by constructing a transfer matrix satisfying the Yang-Baxter equation 27| . 
This matrix was also used in the latter reference to derive the full Yangian symmetry 
algebra of this model, which is ultimately responsible for the highly degenerate character 
of its spectrum. Moreover, the representation theory of the Yangian is closely related to 
Haldane's elegant description of the spectrum in terms of motifs [H3|. It is natural to 
inquire whether a similar construction is possible for the Dn chain of HS type studied in 
this paper. In fact, our numerical calculations show that the spectrum of the Dn chain 
is also highly degenerate, which points to the existence of a large underlying symmetry 
algebra. The characterization of this algebra, and its precise connection with the Yang- 
Baxter equation, is yet another open problem motivated by the present work. 
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